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I.  INTRODUCTION 

This  project  is  concerned  with  the  viscoelastic  properties  of  the 
intervertebral  joint.  Our  previous  studies  under  this  Grant  have  quanti¬ 
fied  the  response  of  the  joint  to  shearing  deformations.  The  work  during 
this  reporting  period  was  directed  toward  the  viscoelastic  response  of 
the  joint  to  torsional  deformations. 

Knowledge  of  the  response  of  the  intervertebral  joints  to  external - 

loading  is  Important  if  the  forces  and  bending  moments  developed  In  the 
spine  are  to  be  related  quantitatively  to  possible  modes  of  injury  or 
abnormal  motions  of  the  joint.  The  use  of  biodynamic  models  has  gained 
widespread  acceptance  as  a  means  of  evaluating  human  tolerance  limits 
and  the  probability  of  injury  In  environments  of  severe  mechanical  stress. 
Numerous  Investigations  have  studied  the  response  of  the  spinal  column 
and  its  constituent  elements  under  axial  tension  and  compression,  in 
torsion,  and  in  bending;  however,  only  a  limited  number  of  investigators 
have  considered  the  response  of  the  spine  to  shear  stress.  The  empirical 
model  of  the  response  of  the  disc  to  a  step  input,  developed  by  Kelley, 


Lafferty  and  Bowman  [1],  is  the  only  attempt  to  provide  an  analytical  repre¬ 
sentation  of  disc  shear  properties.  The  present  study  utilizes  the  same 
experimental  and  modeling  techniques  to  investigate  the  response  of  the 
intervertebral  joint  to  torsional  stress. 

II.  MODEL  ANALYSIS 

Analysis  of  the  response  of  the  intervertebral  to  shear  stress  resulted 
in  an  analytical  representation  of  the  form 

e/o  *  Ej  +  Egt  +  E3  [(1  -  exp  -  tj/tj)  +  (1  -  exp  -  t2/T2)]  (1) 

where  e  =  strain  =  deformation/disc  height 
o  =  stress  =  force/disc  area 
t  *  time 

Ei*  E2,  E3.  tj  and  t2  are  constants 

The  analysis  of  the  shear  response  of  an  elastic  isotropic  material  defines 
the  shear  modulus  G  as 

e/o  *  ^  (2) 

The  analysis  of  the  same  material  subjected  to  torsional  stress  shows  that 

JJ  =  1  (3) 

where  eJ/TH  is  the  torsional  equivalent  e/o  and 

e  ■  angular  deformation  In  radians 
J  =  polar  moment  of  inertia  of  cross  sectional  area  (m*) 

T  =  applied  torque  (N*m) 

H  =  specimen  height  normal  to  the  plane  of  the  applied  moment  (m) 

The  shear  modulus  G  Is  the  same  as  defined  in  Eq.  2. 

Although  the  intervertebral  joint  is  a  complicated  structure  that  Is 
neither  purely  elastic  nor  Isotropic*  some  degree  of  correspondence  between 
shear  response  and  torsional  response  (as  indicated  by  Equations  2  and  3)  is 
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expected  Inasmuch  as  the  disc  structure  responds  to  all  types  of  stress  In 
the  same  basic  mode.  It  was*  therefore*  anticipated  that  the  same  Kelvin 
unit  model  (represented  by  Eq.  1)  developed  for  the  shear  response  would 
characterize  the  torsional  response.  One  would  expect*  however*  that  the 
change  In  force  components  relative  to  the  applied  stress  would  result  In 
different  values  for  the  constants. 

III.  METHODS  AND  MATERIALS 

Test  specimens,  consisting  of  two  vertebrae  and  the  intervening  disc 
with  facets  intact,  were  taken  from  fresh-frozen  Rhesus  monkey  spines.  The 
vertebrae  were  clamped  in  holding  fixtures  with  the  superior  vertebra  con¬ 
strained  to  prevent  angular  displacement  about  the  spinal  axis  but  free  to 
move  in  the  longitudinal  direction.  The  torque  was  applied  to  the  inferior 
vertebra  by  a  MTS  Servoram  controlled  by  a  Cromemco  microprocessor.  The 
torque  was  applied  as  a  step  function  in  increments  of  1.875  Newton-meters 
up  to  failure.  The  creep  response  was  observed  for  90  seconds  at  each  value 
of  torque  and  the  specimen  was  allowed  to  relax  for  a  minimum  of  180  seconds 
between  tests.  One  series  of  tests  (8  specimens)  was  conducted  with  an  axial 
load  of  12.945N  and  a  second  series  (6  specimens)  employed  an  axial  load  of 
27.66N.  The  geometrical  data  (disc  height,  disc  area  and  volume)  were  obtained 
by  sectioning  each  specimen  in  the  longitudinal  direction  and  using  a  Talos 
digitizer  to  record  the  coordinates  of  the  disc.  The  procedures  were  the 
same  as  those  reported  for  the  shear  studies. 

IV.  RESULTS  AND  DISCUSSION 

Typical  responses  of  the  intervertebral  joint  subjected  to  a  constant 
torsional  stress  is  shown  in  Fig.  1  where  the  angular  displacement  (in 


radians)  Is  plotted  vs.  time.  The  plots  show  the  effect  of  torque  level  as 
well  as  that  of  axial  compressive  load.  The  effect  of  Increased  axial  com- 
presslon  Is  more  pronounced  at  the  higher  torque  levels  Indicating  that  the 
articular  surfaces  of  the  facet  joints  are  forced  into  closer  proximity, 
thereby  increasing  resistance  to  increased  angular  deformation. 

Normalization  of  the  data  according  to  Eq.  3  should  permit  the  results 
for  a  given  axial  load  to  be  represented  by  a  function  independent  of  the 
applied  stress  and  specimen  size.  1.e.» 


fg  =  f(kr  k2  —  kr  C,,  C2  —  Cv  t) 

th  2 

where  k^  *  shear  modulus  for  the  il  elastic  element  (N/m  ) 

*  damping  coefficient  for  the  itn  viscous  element  (N*S/m  ) 
t  =  time  (s) 

The  average  values  of  eJ/TH  displayed  a  standard  deviation  on  the  order  of 
±32%  of  the  mean.  We  concluded  that  the  imposition  of  a  specific  axial  load 
to  specimens  obtained  from  subjects  of  widely  varying  size  and  weight  contri¬ 
buted  to  the  magnitude  of  the  standard  deviation.  The  creep  modulus  was, 
therefore,  divided  by  the  axial  load  expressed  as  a  percent  of  body  weight 
of  the  donor.  This  procedure  reduced  the  maximum  standard  deviation  to  the 
order  of  22%.  The  creep  function  is  then  defined  as 


eJ 

THW 


(4) 


where  W  *  (axial  load/body  weight)  x  100 

Use  of  the  model  in  the  form  of  Eq.  1  provided  a  representation  of 
the  results  to  within  ±2%  of  the  averaged  data.  These  results  are  tabulated 
In  Table  I  and  presented  graphically  in  Fig.  2  for  the  two  levels  of  axial 
preload.  The  analytical  model  (represented  by  the  continuous  curve  in 
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Fig.  2)  Is 


K  s  Irw  a  r  *  r +  r [1 ' exp  (‘  +  r  n  -  exP  (-  k2t/c2)i  (5) 

3  3  1  2 

for  which  the  values  of  k  and  C  are  shown  In  Table  II  for  the  torsional  tests 
as  well  as  the  shear  tests. 

Comparison  of  the  model  parameters  Indicate  a  high  degree  of  consistency 
between  the  shear  and  torsional  results ,  l.e. ,  the  time  constants  k/C  are 
identical  In  the  two  cases  and  in  both  cases  kj  *  k2.  These  results  are  not 
surprising  when  one  considers  the  fact  that  the  intervertebral  disc  resists 
all  stress  modes  through  tensile  loading  of  the  fibers  of  the  annulus  fibrosis. 
The  differences  in  the  magnitudes  of  the  corresponding  k  and  C  reflect  (1)  the 
changes  in  the  force  component  directed  along  the  annulus  fibers  as  the  stress 
mode  is  changed  and  (2)  the  effect  of  the  articular  facets  in  the  torsional 
loading. 

The  results  obtained  under  this  grant  provide  an  analytical  model  suitable 
for  use  in  computer  simulations  of  spinal  response  to  shear  and  torsional 
stresses.  The  analytical  representation  Is  not,  however,  unique  in  that 
other  mathematical  functions  (and  mechanical  models)  can  also  represent  the 
temporal  response  of  the  system.  The  model  used  is  not  considered  a  good 
representation  of  the  system  because  increased  response  time  requires  an 
increased  number  of  Kelvin  units  to  accurately  represent  the  response. 

This  results  in  a  series  of  k  and  C  values  to  represent  the  physical  char¬ 
acteristics  of  the  material.  We  are  of  the  opinion  that  a  more  simple 
physical  model  with  time  dependent  elastic  and/or  viscous  coefficients 
exists.  Such  a  model  would  be  more  representative  of  the  actual  material 
and  would  provide  discrete  elastic  and  viscous  element  characterization. 

Our  initial  efforts  to  formulate  such  a  model  are  promising  and  Its  appli¬ 
cation  to  the  general  response  of  the  system  will  be  a  major  part  of  our 
future  investigations. 
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TIME  (seconds 


TABLE  I 

AVERAGE  CREEP  FUNCTION 


12.945  N  Axial  Load  ±  S.E. 

0.8096  ±  0.0319 

I 

0.8508 

0.0332 

0.8728 

0.0346 

0.8958 

0.0351 

0.9105 

0.0359 

0.9636 

0.0383 

1.0161 

0.0401 

1.0512 

0.0429 

1.0908 

0.0422 

1.1174 

0.0418 

1.1376 

0.0444 

1.1516 

0.0458 

1.1675 

0.0461 

1.1742 

0.0471 

Creep  Function,  K(cnr/N' 


27.655  N  Axial  Load  ±  S.E 


26 


12  945  N  Axial  Preload 
27655  N  Axial  Preload 


Fig.  1.  Torsional  Creep  Response  as  a  Function 
of  Applied  Torque,  Axial  Load  and  Time. 


Fig.  2.  Torsional  Creep  Function.  A  Comparison 
of  Model  Response  to  the  Experimental 
Results. 


TABLE  II.  COMPARISON  OF  MODEL  PARAMETERS 
FOR  TORSION  AND  SHEAR  CREEP  TESTS 


Parameter 

Torsion 

Shear 

12.95  N 

Axial  Load 

27.66  N 
Axial  Load 

k1 

5.714  N/mm2 

11.905 

17.39 

k2 

5.714 

11.905 

17.39 

k3 

1.529 

2.865 

1.95 

C1 

4.395  N-S/mn2 

9.158 

13.38 

C2 

43.95 

91.58 

133.8 

C3 

520.92 

795.54 

1000.0 

kl/Cl 

1.3  s'1 

1.3 

1.3 

0.13 


0.13 


0.13 
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